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Abstract. The thermalized energy from the radioactive decays of ^ Ni and Ni and 
their daughter nuclides power the light curves of supernovae near maximum light. The 
bolometric light curve gives us a fundamental understanding of the energy evolution of 
a supernova explosion and the amount of radioactive nuclides produced. In this review, 
I will discuss the bolometric evolution of the Type Hp supernovae SN1987A, and the 
general class of bolometric light curves of Type la thermonuclear explosions. 



1 SN1987A 

We have been monitoring the photometric properties of the Type Hp SN1987A 
in the LMC for more than 15 years. Over the first 5 years the optical photometry 
could be measured from the ground in the typical 1" seeing at our observatories. 
In the near infrared bands of JHK where the seeing is better and the crowding 
stars (which are typically bluer than the SN) are less of a problem, we could 
monitor the SN for up to 10 years. Once the supernova faded below the integrated 
magnitudes of the inner ring which lies at about 1" from the supernova debris, 
HST or ground-based AO imagery is needed to isolate the debris and inner ring 
evolution. In Figure |], I show the structure of the inner ring region. 




Fig. 1. An HST image of region near the inner ring of SN1987A. Region 1 isolates the 
expanding debris. Region 2 represents the shocked region between the debris and the 
inner ring. Region 3 represents the inner ring 
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From our HST images, we have measured the evolution of the brightness 
of the three regions shown in Fig. |l| with respect to an averaged background 
region outside of Region 3. For the earlier ground-based data where we can only 
measure the light from all three regions, we have subtracted off the extrapo- 
lated flux from Regions 2 and 3 to derive estimated fluxes for just the debris 
region. In Fig. ^ we plot the optical light curve for 15 years of evolution. With 
these corrections, the optical and HK photometry from the ground and HST 
flight system equivalents merge well together. The near-infrared photometry in 
J/F1\QW shows a large discontinuity presumably due to the radical differences 
in the two filter systems. The photometry of non-stellar SEDs can show large 
systematic errors due to the differences in the sensitivity functions of the atmo- 
sphere/telescope/filter/detector system. In early photometry of SN1987A, we 
noted difference of up to 0.4 mag in / photometry due to the differences in the 
facility / filters 0. In Table [|, I list the latest photometry for SN1987A 



Table 1. Optical Photometry of SN1987A as of May 2002 





U 


B 


V 


R 


/ 


debris 


20.91(32) 


20.88(12) 


21.28(07) 


20.96(22) 


20.46(04) 


decay rate (mag y^^) 


0.23 


0.20 


0.24 


0.35 


0.22 



Table 2. Near-Infrared Photometry of SN1987A 



days since explosion 


J 


H 


K 


4151 


19.81(09) 


18.29(06) 


18.40(02) 


5428 


20.46(15) 







The leveling off of the light curves after year 5 is due to two causes: the 
longer decay times for the remaining radioactive nuclides and the freeze-out of 
the cooling. By year five, the main radioactive energy sources of ^^Ni and its 
daughter nucleus ^^Co with e-folding times of 8.8d and 111.3d have decayed 
away. The nuclides ^^Co (the daughter nuclide of ^^Ni) and ''^Ti with e-folding 
times of 390d and 87y are left as the main energy input, along with the energy 
input from a possible (but as yet unseen) pulsar. However, the light curves at 
this phase no longer represent the prompt thermalization of the input radioactive 
energy. As shown in ||J^, the recombination and cooling time scales become 
larger than the expansion time scale, and the debris nebula is not able to cool 
at the rate of radioactive energy input. 

The ultraviolet, optical, and infrared photometry can be integrated to the 
thermalized "uvoir" energy flux to derive the masses of the synthesized radioac- 
tive nuclides. The early time data (less than lOOOd) with prompt thermalization 



Bolometric Light Curves of Supernovae 3 




a e> ID ifi 5 10 15 

Years years 

Fig. 2. Optical light curves (left panel) and near-infrared light curves (right panel) 
for the debris (Region 1) in SN1987A. The data have been shifted by the magnitudes 
listed in the legend. The data are a combination of ground-based UBVRIJHK and 
the equivalent HST flight system magnitudes. 



and cooling can be simply fit to radioactive decay models as in . The late- 

time bolometric light curves, when fit with models including the freeze-out, can 
be used to measure the amount of '"^^Ni and "^^Ti. Fransson and Kozma (|^) 
present the latest values for the amounts of these nuclides where they found 
the following masses: M(57Ni,5^Ni,'''*Ti)=(O.O69,O.OO3,O.OOOl)M0. Lundqvist et 
al. (0]) find upper limits on the mass of ^■^Ti consistent with these numbers 
based on the non-detection of the 25m/x lines of [Fel] and [Fell] from ISO/SWS. 
They caution, however, that if dust cooling is important, the limits on the mass 
of '*''Ti could be significantly higher. Only direct detection of the 1.157 MeV line 
of 4'^Ti will resolve the ambiguity in the models. 

In Fig. ^ I show the late-time photometric data for the ring and the inter-ring 
region. The ring region suddenly began to brighten at a rate of ~ —0.24 mag 
y~^ in UBVI and ~ —0.12 mag y^^ in R since year 13. The brightening of the 
ring was predicted in The ring was expected to brighten by 7.5mag to 

about = 11 as the blast wave from the debris struck the inner ring starting 
around year 16-20. However, in year 10, Pun et al. ((ll)) discovered a single 
spot brightening, evidently due to an inward protrusion of the inner ring. Our 
observations here show that the general rapid brightening of the ring began later, 
around year 13. 
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The inter- ring region has accelerated its brightening since year 13 and is 
now brightening at a rate of ~ —0.12 mag in UBVI. The inter-ring light 
is presumably due to the the emission from the reverse shock which is located 
at ~ 75% of the radius of the inner boundary of the inner ring. This emission 
was predicted by ||l^ and discovered by ||l^ with STIS data from HST. Refined 
models are presented in [ p^p^ . The reverse shock formed in the equatorial 
plane as neutral hydrogen atoms streaming from the debris at 15000 km s~^ 
hit the ionized region at the inner surface of the ring. There is also diffuse light 
present in the STIS data which may be due to excitation by non-thermal particles 
accelerated by the shock. 




Fig. 3. UBVRI optical light curves for the ring Region 3 (left panel) and the inter-ring 
Region 2 (right panel). The magnitudes from bottom to top are BVIUR. 



Significant dust formed in SN1987A starting around day 500, and by day 
1000, the bolometric flux of the nebula was dominated by dust radiation with 
a temperature of Tdust ~ 200K and mid-infrared line emission(fl|l5t). In Fig. || 

1 show the latest data on the mid-infrared evolution of SN1987A with the ISO- 
CAM detection (Q) and a detection from OSCIR on the CTIO 4m telescope 
([T^) at 10|J,m. Without 20|a.m detections however, we can't estimate the dust 
temperature and cannot measure the dust emission. If we assume the same 
bolometric corrections from day 1800, we find the bolometric flux of SN1987A 
is logio(L) ~ 36.0 erg s~^. 

2 The Bolometric Behavior of Type la Supernovae 

Because there is much less envelope mass above the radioactive energy sources 
in Type la SNe, the bolometric properties of these supernovae are more difficult 
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Fig. 4. Mid-infrared light curves of SN1987A 



to interpret because much less of the radioactive energy is thermalized. The 
physics of the thermalization is discussed in . The fraction of the radiation 

which is thermalized must be transported to the surface and radiated away. If 
the diffusion time for the energy is short compared to the dynamical time, the 
energy appears as thermalized "uvoir" radiation. If the diffusion time is long, the 
trapped radiation will be converted partly into kinetic energy and not seen in the 
bolometric light curve. With a longer diffusion time, the peak of the bolometric 
light curve is shifted to later times when the input energy from the radioactive 
decays is less. This leads to "Arnett's Law" which states that the bolometric 
luminosity at maximum light is equal to the instantaneous energy released by 
the radioactive decay. This law was (and is) an important tool to measure the 
amount of radioactive nickel produced in the explosion. 

Continuing the discussion in , for the first 60 days since explosion the ther- 
malization time is short compared to the dynamical time, and little radioactive 
energy is converted into kinetic energy. Thus the integrated uvoir bolometric 
luminosity represents the fraction of radioactive energy input which has been 
thermalized. By day 40 more than one half of gamma-ray radiation is leaving 
the explosion unthermalized. Even at B maximum light some 20 days after ex- 
plosion, some 15% of the gamma rays are unthermalized. The observed uvoir 
bolometric light curve is then a function of the rapidly decreasing optical depth 
to the gamma-rays over the first 100 days and the decreasing energy input from 
the radioactive nickel nuclides. The combination of these two processes acceler- 
ates the light curve decay at a rate faster than the radioactive decay of ^^Co 
after maximum light. 

The uvoir bolometric light curve interpretation is also complicated by the 
efficiency of the trapping of the gamma rays and positrons in the ^^Co decay 
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(|2^,^,|§|). Milne et al. ((2l)) show that the gamma-rays from the '^''Co decay 
carry 30 times more energy than the positrons, but also are more penetrating 
and can escape more easily. As the nebula expands from days 50 to 200, the 
energy deposition switches from the '^^Co gamma-ray Comptonization to the 
energy input from the positron annihilations. Some of the positrons, however, 
can be transported out of the nebula, depending on the magnetic fields structure. 
The bolometric light curve at this phase then becomes dependent on the details 
of the mixing of the ^^Co in the nebula and the properties of the magnetic fields. 

While the interpretation of the bolometric light curve requires theoretical 
knowledge of the optical depths to gamma-rays and the details about the mix- 
ing of ^^Ni, the construction of the uvoir bolometric light curve is rather simple 
due to the following coincidence - - most of the thermalized flux appears at op- 
tical wavelengths. In Fig. ^ I show the cumulative fiux of SN1992A at various 
phases during the first 100 days of evolution ([§4j). Also shown in the figure is a 
comparison of a Type la SED with that of a Type Ic and II ([^,^). It can be 
clearly seen that the broad peak of the flux distribution for Type la SNe appears 
in the optical region. In fact, typically 80% or more of the uvoir flux appears in 
the optical from day -6 onwards. 

Only a few papers have been published trying to estimate the bolometric light 
curves (see for a summary). In ]27| , p8t , the optical broad-band magnitudes 
were integrated to provide a magnitude that should be a close surrogate to the 
uvoir bolometric magnitude. From the estimated bolometric light curve, they 
found a range of more than a factor in 10 in the ^^Ni masses for a group of nearby 
Type la SNe. In j2^, a V magnitude was used with a bolometric correction to 
study the gamma-ray trapping in the late-time light curves, which also showed 
a significant range in ^^Ni masses. 

In 1996, I used the small amount of data on supernovae which had space 
ultraviolet spectra, optical UBVRI, and near-infrared JHK data to estimate 
accurate bolometric fluxes (||2^). I found that there was a range in peak bolo- 
metric magnitudes (implying a range in nickel masses) and that the bolometric 
light curves appear to have a small secondary hump in the light curve around 
days 20-40 corresponding to the secondary maximum in the / band. Such a 
flux redistribution was unexpected, and points to a rather sudden change in the 
opacity and cooling in the nebula. 

The construction of the bolometric light curve requires the sum of the space 
ultraviolet, optical, and near-infrared data. The falloff of the flux beyond H is 
such that the assumption of a R-J law adds only a few percent to the total flux 
and a R-J extrapolation from H or K is entirely adequate. The extrapolation 
to the space ultraviolet, however, is larger. Using 1992A data, I found that 
extrapolating the U flux point at 360nm to zero flux at 300nm was a reasonable 
representation of the true space ultraviolet flux from day -6 onwards, yielding 
agreement to 5% in the uvoir flux between the full integration and the integration 
using the ultraviolet extrapolation. It is important to continue to observe nearby 
SNe with HST in the ultraviolet to understand the diversity of the ultraviolet 
flux, which may be an indicator of the metallicity of the progenitor. 
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Fig. 5. Flux distributions of Type la SNe. The left panel shows the cumulative flux 
distribution for SNf992A for days -6,0,5,20,80 (from top to bottom at A = 300nm), 
with the distribution for day shown in the darker line. 80% or more of the total uvoir 
flux appears in the optical window of 300-lOOOnm. The right panel shows the SEDs 
(top to bottom) for SN1999em (Type II), SN 1999ex (Ibc) and SN 1999ee (Type la) 
near maximum light (reference Pq]). Note that the broad flux peak of the SED appears 
in the optical. 

Our groups at CTIO and LCO have been following nearby Type la and Type 
II in UBVRIYJHK in part to study the bolometric properties of these super- 
novae. The new filter Y at l)im ([^) has been designed to fit in a wavelength 
band which is remarkably free of telluric absorption and provides a flux point 
between the widely separated / and J filters. Figure || ([Q) shows the SED 
evolution of the nearby SN 2001el. This figure shows that the SED is well fit 
by a black-body with T^ff = 14,000K before maximum light steepening to a 
R-J law by maximum light in BVRIJHK. After maximum light, a flux deflcit 
appears in J and the the bands U BV drop away from a thermal distribution. 

In Figure |^ I plot the uvoir bolometric light curves calculated from the inte- 
gration of the broad-band magnitudes, including an extrapolation to the ultra- 
violet to account for the unobserved space ultraviolet, and to the infrared using 
a R-J law extending from / or K. The absolute flux scale has been set by using 
either SBF, Cepheid, TRGB, or PN distances for the nearby SNe, or a Hub- 
ble law of Hq — 74 km s^^ Mpc~^ for the more distant SNe. The comparison 
of SN1999ee and SN2001el shows that the location and size of the secondary 
bump is significantly different despite the fact that the peak luminosities are 
very similar. 

Also plotted in Fig. are the uvoir bolometric light curves for 16 Type la 
SNe. All have UBVRI photometry and 10 have JHK photometry. The sense of 
this figure is that most SNe have roughly the same peak bolometric luminosity, 
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Fig. 6. SED evolution for the Type la SN2001el. The left panels shows the SEDs for 
days (since Bmax) -10.7 and +1.3. The right panel shows the SEDs for days 1.3, 13.3, 
25.2, 40.2, and 64.3. ATeff = 14,000K blackbody is fit to the -10.7 and 64.3 data, 
and a Rayleigh- Jeans law is fit to the data 1.3 data. The SED steepens to the blue as 
the SN approaches maximum light and then reddens after maximum. The fiux deficit 
in J appears right after maximum light. Note that the IHK data roughly follow the 
R-J law. 

but the size and placement of the secondary hump is quite variable. For instance, 
the peculiar SN2000cx has a rather normal bolometric light curve near maximum 
light, but almost completely lacks a secondary hump at 30 days. This bolomet- 
ric behavior reflects the fact that SN2000cx had an anomalously weak /-band 
secondary maximum. Evidently the peak bolometric magnitude (and thus the 
^^Ni mass) is not strongly coupled to the opacity and flux redistribution causing 
the secondary hump. 

Finally, in Fig. |^ I show the peak bolometric luminosity versus the light curve 
shape parameter Ami5 for the sample plotted in Fig. |^. Rather than the roughly 
linear relationship between the intrinsic luminosity in i? or 1/ and Amis, we find 
that for much of the range of Amis there is little or no relationship between 
logio(L) at peak and Amis, at least for Amis < 1-3. 
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Fig. 7. "UVOIR" bolometric light curves for Type la SNe. These light curves include 
extrapolations to the space ultraviolet and the extrapolation to the red using a R- 
J law. Closed circles indicate an extrapolation from K and open circles indicate an 
extrapolation from /. The left panel shows the bolometric light curves for SN2001el 
(upper curve, Amis = 1-13) and SN1999ee (lower curve, Amis = 0.94). The SN1999ee 
data have been shifted by 0.25dex for presentation purposes. The right panel shows 
the bolometric light curves for 16 Type la SNe. 
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Fig. 8. The pea k uvoir bolometric luminosity plotted against the light curve shape 
parameter ol for the SNe shown in Fig. |^. 
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